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□ 1. Document ID: US 20030049867 Al 
Using default format because multiple data bases are involved. 

L16: Entry 1 of 11 File: PGPB 



Mar 13, 2003 



PGPUB- DOCUMENT-NUMBER: 20030049867 
PGPUB- FILING-TYPE: new 

DOCUMENT-IDENTIFIER: US 20030049867 Al 

TITLE: Methods for parallel detection of compositions having desired 

characteristics 

PUBLICATION-DATE: March 13, 2003 



INVENTOR-INFORMATION: 

NAME 

Selifonov, Sergey A. 
Huisman, Gjalt W. 



CITY 

Plymouth 
San Carlos 



STATE 

MN 

CA 



COUNTRY 

US 

US 



US-CL-CURRENT: 436 / 518 ; 435 /6, 702 /19 




□ 2. Document ID: US 20030022105 Al 

L16: Entry 2 of 11 File: PGPB 



Jan 30, 2003 



PGPUB- DOCUMENT -NUMBER : 20030022105 
PGPUB-FILING-TYPE: new 

DOCUMENT-IDENTIFIER: US 20030022105 Al 

TITLE: TWO -PHOTON UPCONVERTING DYES AND APPLICATIONS 
PUBLICATION-DATE: January 30, 2003 



I NVENTOR- I N FORMAT I ON : 
NAME 

PRASAD, PARAS N. 
BHAWALKAR, JAY ANT D. 
CHENG, PING 
PAN, SHAN JEN 



CITY 
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□ 3. Document ID: US 20010024796 A 1 

L16: Entry 3 of 11 File: PGPB 



Sep 27, 2001 



PGPUB-DOCUMENT-NUMBER: 20010024796 

PGPUB- FILING-TYPE : new 

DOCUMENT-IDENTIFIER; US 20010024796 Al 

TITLE: Methods for parallel detection of compositions having desired 
characteristics 

PUBLICATION-DATE: September 27, 2001 



INVENTOR-INFORMATION : 
NAME 

Selifonov, Sergey A. 
Huisman, Gjalt W. 



CITY 

Gainesville 
San Carlos 



STATE 

FL 

CA 



COUNTRY 

US 

US 



US-CL-CURRENT: 435/7.1/ 324/308 




□ 4. Document ID: US 6891371 Bl 
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File: USPT 



May 10, 2005 



US-PAT-NO: 6891371 

DOCUMENT-IDENTIFIER: US 6891371 Bl 

** See image for Certificate of Correction ** 

TITLE: Method and system of generating an MRS spectrum from multiple receiver data 
DATE-ISSUED: May 10, 2005 



INVENTOR- INFORMATION: 






NAME 


CITY 


STATE 


Frigo; Frederick J. 


Waukesha 
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Heinen; James A. 


Wauwatosa 
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Raidy; Thomas E. 


Elm Grove 
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Hopkins; Jeffery A. 


Pewaukee 
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ZIP CODE 
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□ 5. Document ED: US 6402037 Bl 

L16: Entry 5 of 11 



File: USPT 



Jun 11, 2002 



US-PAT-NO: 6402037 

DOCUMENT-IDENTIFIER: US 6402037 Bl 

TITLE: Two-photon upconverting dyes and applications 
DATE-ISSUED: June 11, 2002 



INVENTOR-INFORMATION: 
NAME 

Prasad; Paras N. 
Bhawalker; Jayant D. 
Cheng; Ping Chin 
Pan; Shan Jen 



CITY 


STATE 


Williamsville 
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□ 6. Document ID: US 6037772 A 
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File: USPT 



Mar 14, 2000 



US- PAT-NO: 6037 772 

DOCUMENT-IDENTIFIER: US 6037772 A 
TITLE: Fast spectroscopic imaging system 
DATE-ISSUED: March 14, 2000 



INVENTOR-INFORMATION: 
NAME 

Karczmar; Gregory S. 
Kovar; David A. 



CITY 

Crete 

Chicago 



STATE 

IL 

IL 



ZIP CODE 



COUNTRY 



US-CL-CURRENT: 324/309; 324/318, 600/410 



□ 7. Document ED: US 5912257 A 

L16: Entry 7 of 11 



File: USPT 



Jun 15, 1999 



US-PAT-NO: 5912257 

DOCUMENT- IDENTIFIER: US 5912257 A 



TITLE: Two-photon upconverting dyes and applications 
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DATE-ISSUED: June 15, 1999 



I^^VENTOR-INFORMATION : 
NAME 

Prasad; Paras N. 
Bhawalkar; Jayant D. 
He; Guang S. 
Zhao; Chan F. 
Gvishi; Raz 
Ruland; Gary E. 
Zieba; Jaroslaw 
Cheng; Ping Chin 
Pan; Shan Jen 



CITY 
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Tonawanda 
Williamsville 
San Diego 
K. Tiron 
Grand Island 
Santa Rosa 
Williamsville 
Amherst 



STATE 

NY 

NY 

NY 

CA 

NY 
CA 
NY 
NY 



ZIP CODE 



COUNTRY 



IL 



US -CL- CURRENT: 514/356; 250 / 338 . 1 , 430/338, 430/343, 514/709, 522/6, 546/329, 
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□ 8. Document ID: US 5903 149 A 

L16: Entry 8 of 11 



File: USPT 



May 11, 1999 



US-PAT-NO: 5903149 

DOCUMENT-IDENTIFIER: US 5903149 A 

** See image for Certificate of Correction ** 

TITLE: Three-dimensional localized proton NMR spectroscopy using a hybrid of one- 
dimensional hadamard with two-dimensional chemical shift imaging 

DATE-ISSUED: May 11, 1999 



INVENTOR-INFORMATION : 
NAME 

Gonen; Oded 
Goelman; Gadi 
Leigh; John S. 
Bolinger; Lizann 



CITY 

Cheltenham 
Rosh Haaein 
Philadelphia 
Philadelphia 



STATE 
PA 

PA 
PA 



ZIP CODE 



COUNTRY 



IL 



US-CL-CURRENT : 324/307 





□ 9. Document ID: US 5570019 A 

L16: Entry 9 of 11 



File: USPT 



Oct 29, 1996 



US-PAT-NO: 5570019 

DOCUMENT- IDENTIFIER: US 5570019 A 
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TITLE: Method for magnetic resonance spectroscopic imaging with multiple spin- 
echoes 

DATE-ISSUED: October 29, 1996 



INVENTOR-INFORMATION: 
NAME 

Moonen; Chrit T. W. 
Duyn; Jeff 



CITY 

Silver Spring 
Kensington 



STATE 

MD 

MD 



ZIP CODE 



COUNTRY 



US-CL-CURRENT: 324/309; 324/307 



□ 10. DocumentID: US 5192909 A 

L16: Entry 10 of 11 



File: USPT 



Mar 9, 1993 



US- PAT-NO: 5192909 

DOCUMENT-IDENTIFIER: US 5192909 A 

TITLE: Spectroscopic localization using pinwheel NMR excitation pulses 
DATE-ISSUED: March 9, 1993 



INVENTOR-INFORMATION: 
NAT^E 

Hardy; Christopher J. 
Bottomley; Paul A. 
Cline; Harvey E. 



CITY 

Schenectady 
Clifton Park 
Schenectady 



STATE 
NY 
NY 
NY 



ZIP CODE 



COUNTRY 



US-CL-CURRENT: 324/309; 324/307 



□ 11. DocumentID: US 4585992 A 

L16: Entry 11 of 11 File: USPT Apr 29, 1986 

US-PAT-NO: 4585992 

DOCUMENT-IDENTIFIER: US 4585992 A 

TITLE: NMR imaging methods 

DATE-ISSUED: April 29, 1986 

INVENTOR-INFORMATION : 
NAME 

Maudsley; Andrew A. 
Hiial; Sadek K, 



CITY STATE ZIP CODE COUNTRY 

Woburn MA 
New York NY 
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Search Results - Record(s) 1 through 4 of 4 returned. 



□ 1. DocumentID: US 6891371 Bl 
Using default format because multiple data bases are involved. 



L17: Entry 1 of 4 



File: USPT 



May 10, 2005 



US- PAT-NO: 68 91371 

DOCUMENT-IDENTIFIER: US 6891371 Bl 

** See image for Certificate of Correction ** 

TITLE: Method and system of generating an MRS spectrum from multiple receiver data 
DATE-ISSUED: May 10, 2005 



INVENTOR- INFORMATION: 






NAME 


CITY 


STATE 


Frigo; Frederick J. 


Waukesha 


WI 


Heinen; James A. 


Wauwatosa 


WI 


Raidy; Thomas E . 


Elm Grove 


WI 


Hopkins; Jeffery A. 


Pewaukee 


WI 



ZIP CODE 



COUNTRY 



US-CL-CURRENT: 324/307; 324/309 




□ 2. DocumentID: US 5570019 A 

L17: Entry 2 of 4 



File: USPT 



Oct 29, 1996 



US-PAT-NO: 5570019 

DOCUMENT-IDENTIFIER: US 5570019 A 

TITLE: Method for magnetic resonance spectroscopic imaging with multiple spin- 
echoes 

DATE-ISSUED: October 29, 1996 



INVENTOR-INFORMATION : 
NAME 

Moonen; Chrit T. W. 
Duyn; Jeff 



CITY 

Silver Spring 
Kensington 



STATE 

MD 

MD 



ZIP CODE 



COUNTRY 



US-CL-CURRENT: 324/309; 324/307 
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□ 3 . Document ID: US 5 1 92909 A 

L17: Entry 3 of 4 



File: USPT 



Mar 9, 1993 



US-PAT-NO: 5192909 

DOCUMENT-IDENTIFIER: US 5192909 A 

TITLE: Spectroscopic localization using pinwheel NMR excitation pulses 
DATE-ISSUED: March 9, 1993 



INVENTOR-INFORMATION: 
NAME 

Hardy; Christopher J. 
Bottomley; Paul A. 
Cline; Harvey E. 



CITY 

Schenectady 
Clifton Park 
Schenectady 



STATE 
NY 
NY 
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ZIP CODE 
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US-CL-CURRENT: 324/309; 324/307 



□ 4. Document ID: US 4585992 A 

L17: Entry 4 of 4 



File: USPT 



Apr 29, 1986 



US-PAT-NO: 4585992 

DOCUMENT-IDENTIFIER: US 4585992 A 
TITLE: NMR imaging methods 
DATE-ISSUED: April 29, 1986 



INVENTOR-INFORMATION : 
NAME 

Maudsley; Andrew A. 
Hilal; Sadek K. 
Simon; Howard E. 
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Woburn 
New York 
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L17: Entry 1 of 4 



File: USPT 



May 10, 2005 



DOCUMENT-IDENTIFIER: US 6891371 Bl 

** See image for Certificate of Correction ** 

TITLE: Method and system of generating an MRS spectrum from multiple receiver data 



Abstract Text (1) : 

A system and method for multi -channel MR spectroscopy (MRS) includes simultaneously 
acquiring MR signals from multiple coils and processing the MR signals individually 
to generate multiple sets of MRS results. The multiple sets of MRS results are 
combined to form a single superset of MRS results that is displayed as a simple 
proton MRS absorption spectrum. 

Brief Summary Text ( 2 ) : 

The present invention relates generally to magnetic resonance spectroscopy (MRS) 
and, more particularly, to a system and method for multiple- receiver proton 
spectroscopy such that a single absorption spectrum is generated as a combination 
of data received from multiple receiver coils . 

Brief Summary Text (3) : 

When a substance such as human tissue is subjected to a uniform magnetic field 
(polarizing field B.sub.O), the individual magnetic moments of the spins in the 
tissue attempt to align with this polarizing field, but precess about it in random 
order at their characteristic Larmor frequency. If the substance, or tissue, is 
subjected to a magnetic field ( excitation field B.sub.l) which is in the x-y plane 
and which is near the Larmor frequency, the net aligned moment, or "longitudinal 
magnetization", Mz, may be rotated, or "tipped", into the x-y plane to produce a 
net transverse magnetic moment Mt. A signal is emitted by the excited spins after 
the excitation signal B.sub.l is terminated and this signal may be received and 
processed to form an image or a magnetic resonance spectroscopy . 

Brief Summary Text (4) : 

When utilizing these signals to produce images, magnetic field gradients (G.sub.x, 
G.sub.y, and G.sub.z) are employed. Typically, the region to be imaged is scanned 
by a sequence of measurement cycles in which these gradients vary according to the 
particular localization method being used. The resulting set of received NMR 
signals is digitized and processed to reconstruct the image using one of many well 
known reconstruction techniques. 

Brief Summary Text (5) : 

Magnetic Resonance Spectroscopy (MRS) may be used in vivo for the determination of 
individual chemical compounds located within a volume of interest. The underlying 
principle of MRS is that atomic nuclei are surrounded by a cloud of electrons which 
slightly shield the nucleus from any external magnetic field. As the structure of 
the electron cloud is specific to an individual molecule or compound, the magnitude 
of this screening effect is then also a characteristic of the chemical environment 
of individual nuclei. Since the resonant frequency of the nuclei is proportional to 
the magnetic field it experiences, the resonant frequency can be determined not 
only by the external applied field, but also by the small field shift generated by 
the electron cloud. Detection of this chemical shift, which is usually expressed as 
"parts per million" (PPM) of the main frequency, requires high levels of 
homogeneity of the main magnetic field Be. 
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Brief Summary Text { 6) : 

Typically, MR proton spectroscopy is used to generate a one-dimensional (ID) 
frequency spectrum representing the presence of certain chemical bonds in the 
region of interest. In medical diagnosis and treatment, MRS provides a non-invasive 
means of identifying and quantifying metabolites from a region of interest, often 
the human brain. By finding the relative spectral amplitudes resulting from 
frequency components of different molecules, medical professionals can identify 
metabolites indicative of diseases, disorders, and other pathologies such as 
Alzheimer's disease, cancer, stroke, and the like. In this context, two nuclei are 
typically of particular interest, HR and p. sup. 31. Phosphorus 31 MRS is directed to 
the detection of compounds involved in energy metabolism relating to membrane 
synthesis and degradation. Metabolites of particular interest in proton MRS studies 
include glutamate (Glu) , glutainine (Gin), choline (Cho), creatine (Cre) , N- 
acetylaspartate (NAA) , and the inositols (ml and si). 

Brief Summary Text ( 7 ) : 

Recently, increased sampling resolutions and the ability to simultaneously receive 
data fro m multiple -channel receiver coils has lead to greater accuracy in MR data 
acquisition. Accordingly, radiologists and other medical personnel are able to 
review separate MRS results from each channel for data acquired during a single 
scan. Additionally, improved data acquisition capabilities have created 
opportunities for implementing new signal processing algorithms to improve the 
sensitivity and accuracy of MRS procedures. Specifically, nonparametric techniques, 
such as Capon analysis and Amplitude and Phase Estimation of a Sinusoid (APES) 
analysis can provide 2D results for MRS experiments that show both frequency and 
damping information. These techniques provide improved sensitivity and 
differentiation between metabolites . 

Brief Summary Text (8) : 

However, whether implementing ID or 2D proton spectroscopy, interpreting the 
results of multiple channels can be a tedious and time consuming process because a 
medical professional must inspect the results from each channel independently. That 
is, since a separate result for each channel is generated, when reviewing the 
results, the medical professional must independently interpret the results from 
each channel. This process of reviewing multiple sets of results and changing 
therebetween creates an arduous review process that is susceptible to erroneous 
human estimations and decreases patient throughput. That is, since the medical 
professional reviewing the results must repeatedly switch between the results from 
each independent channel, the medical professional is required to estimate an 
overall review of data acquired from the various independent channels. For example 
in an 8 channel head coil, a medical professional must individually review data 
from 8 different receiver channels to make a diagnosis. Relying on human estimation 
not only lengthens the review process but also complicates the diagnostic process. 

Brief Summary Text (9) : 

It would therefore be desirable to have a system and method capable of generating a 
single set of results from a multiple receiver MR proton spectroscopy scan. 

Brief Summary Text (11) : 

The present invention provides a system and method of generating a single 
compilation of results from an MR proton spectroscopy scan that overcomes the 
aforementioned drawbacks. Specifically, the present invention provides a system and 
method of processing MR signals individually and then combining the multiple sets 
of MRS results to generate and display a single compilation of MRS results. As 
such, a radiologist or other medical personnel may make a diagnosis from a single 
composite spectrum . 

Brief Summary Text (12): 

In accordance with one aspect of the invention, a method of multi -channel MR 
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spectroscopy (MRS) is disclosed that includes simultaneously acquiring MR signals 
from multiple coils and processing the MR signals individually to generate multiple 
sets of MRS results. The method also includes combining the multiple sets of MRS 
results to form a single superset of MRS results and generating and displaying a 
proton MRS absorption spectrum from the single superset of MRS results for clinical 
inspection . 

Brief Summary Text (13) : 

In accordance with another aspect of the invention, a magnetic resonance imaging 
(MRI) apparatus includes a plurality of gradient coils positioned about a bore of a 
magnet to impress a polarizing magnetic field and is used to generate a set of MRS 
results. An RF transceiver system and an RF switch are controlled by a pulse module 
to transmit and receive RF signals to and from an^ RF coi l assembly to acquire MR 
images. The MRI apparatus also includes a compufeT*^ programmed to acquire MR signals 
from multiple receiver channels and process the MR signals from each channel 
independently. The computer is further programmed to process MR signals from each 
coil into a single set of MRS results for each coil^ then combine these results 
into a single composite set of MRS results and display a proton MRS absorption 
spectrum from the single composite set of MRS results. 

Brief Summary Text (14) : 

In accordance with another aspect of the invention, the invention is embodied in a 
computer program stored on a computer readable storage medium and having 
instructions which, when executed by a computer, cause the computer to process MRS 
results from the multiple coils independently. The computer combines the processed 
MRS results into a single composite set of MRS results and displays the set as a 
single MRS absorption spectrum . 

Drawing Description Text (5) : 

FIG. 2 is a flow chart setting forth the steps of compiling a single absorption 
spectrum from a multiple receiver MR proton spectroscopy scan in accordance with 
one embodiment of the present invention. 

Drawing Description Text (6) : 

FIG. 3 is an example of a single ID absorption spectrum compiled from a multiple 
receiver MR proton spectroscopy scan in accordance with one embodiment of the 
present invention. 

Drawing Description Text ( 7 ) : 

FIG. 4 is an example of a single 2D absorption spectrum 
receiver MR proton spectroscopy scan in accordance with 
present invention. 

Detailed Description Text (2) : 
Referring to FIG. 1, the major components of a preferred MRI system 10 
incorporating the present invention are shown. The operation of the system is 
controlled from an operator console 12 which includes a keyboard or other input 
device 13, a control panel 14, and a display screen 16. The console 12 communicates 
through a link 18 with a separate computer system 20 that enables an operator to 
control the production and display of images and MRS results on the display screen 
16. The computer system 20 includes a number of modules which communicate with each 
other through a backplane 20a. These include an image processor module 22, a CPU 
module 24 and a memory module 26, known in the art as a frame buffer for storing 
image data arrays . The computer system 20 is linked to disk storage 28 and magneto 
optical disk (MOD) drive 30 for storage of image data and programs, and 
communicates with a separate system control 32 through a high speed serial link 34. 
The input device 13 can include a mouse, joystick, keyboard, track ball, touch- 
activated screen, light wand, voice control, or any similar or equivalent input 
device, and may be used for interactive geometry prescription. 



compiled from a multiple 
one embodiment of the 
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Detailed Description Text (3) : 

The system control 32 includes a set of modules connected together by a backplane 
32a. These include a CPU module 36 and a pulse generator module 38 which connects 
to the operator console 12 through a serial link 40. It is through link 40 that the 
system control 32 receives commands from the operator to indicate the scan sequence 
that is to be performed. The pulse generator module 38 operates the system 
components to carry out the desired scan sequence and produces data which indicates 
the timing, strength and shape of the RF pulses produced, and the timing and length 
of the data acquisition window. The pulse generator module 38 connects to a set of 
gradient amplifiers 42, to indicate the timing and shape of the gradient pulses 
that are produced during the scan. The pulse generator module 38 can also receive 
patient data from a physiological acquisition controller 44 that receives signals 
from a number of different sensors connected to the patient, such as ECG signals 
from electrodes attached to the patient. And finally, the pulse generator module 38 
connects to a scan room interface circuit 46 which receives signals from various 
sensors associated with the condition of the patient and the magnet system. It is 
also through the scan room interface circuit 46 that a patient positioning system 
48 receives commands to move the patient to the desired position for the scan. 

Detailed Description Text (4) : 

The gradient waveforms produced by the pulse generator module 38 are applied to the 
gradient amplifier system 42 having G.sub.x, G.sub.y, and G.sub.z amplifiers. Each 
gradient amplifier excites a corresponding physical gradient coil in a gradient 
coil assembly generally designated 50 to produce the magnetic field gradients used 
for spatially encoding acquired signals. The gradient coil assembly 50 forms part 
of a magnet assembly 52 which includes a polarizing magnet 54 and a whole -body RF 
coil 56. A transceiver module 58 in the system control 32 produces pulses which are 
amplified by an RF amplifier 60 and coupled to the RF coil 56 by a transmit/ receive 
switch 62. The resulting signals emitted by the excited nuclei in the patient may 
be sensed by the same RF coil 56 and coupled through the transmit/ receive switch 62 
to a preamplifier 64. The amplified MR signals are demodulated, filtered, and 
digitized in the receiver section of the transceiver 58. The transmit /receive 
switch 62 is controlled by a signal from the pulse generator module 38 to 
electrically connect the RF amplifier 60 to the coil 56 during the transmit mode 
and to connect the preamplifier 64 to the coil 56 during the receive mode. The 
transmit /receive switch 62 can also enable a separate RF coil {for example, a 
surface coil) to be used in either the transmit or receive mode. 

Detailed Description Text (5) : 

The MR signals picked up by the RF coil 56 are digitized by the transceiver module 
58 and transferred to a memory module 66 in the system control 32. A scan is 
complete when an array of MR raw data has been acquired in the memory module 66. 
This MR raw data is input to an array processor 68 which processes the MR raw data 
as necessary to create the desired results which may include image data or MRS 
results. The image data or MRS results are conveyed through the serial link 34 to 
the computer system 20 where they are stored in memory, such as disk storage 28. In 
response to commands received from the operator console 12, the image data or MRS 
results may be archived in long term storage, such as on the MOD drive 30, or may 
be further processed by the image processor 22 and conveyed to the operator console 
12 and presented on the display 16. 

Detailed Description Text (6) : 

The present invention includes a system and method suitable for use with the above- 
referenced MR system, or any similar or equivalent system for obtaining MR data. 
Proton MRS can be implemented on the above-referenced MR system and is used in vivo 
to measure the concentration of a number of metabolites . High-field magnetic 
resonance scanners (1.5 T or greater) are typically used to perform MRS studies 
involving a number of pathologies, including but not limited to, cancer, multiple 
sclerosis, Alzheimer's disease, and acute stroke. 
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Detailed Description Text (7) : 

Referring now to FIG. 2, an MRS scan is conducted to simultaneously acquire data 
for a ID or 2D absorption spectrum plot from multiple receiver channels of an MR 
system. The MRS scan includes a reference data acquisition in which a signal is 
acquired from a region or volume of interest. The MRS scan also includes a water - 
suppressed data acquisition in which the dominant water signal has been removed via 
a chemical selective saturation, such as a known CHESS pulse sequence. The MRS data 
acquisition requires that a number of scan parameters be determined for each region 
of interest. The center frequency, flip angle, transmit gain, receive gain, and 
shimming parameters are determined for each MRS experiment during prescan . These 
parameters are optimized for a particular TE and TR selected for a particular scan. 

• 

Detailed Description Text (9) : 

where A.sub.r, . theta . . sub . r, . sigma . . sub . r, and . omega sub . r, are the amplitude, 
phase, damping factor and angular frequency of the r.sup.th component, 
respectively, and .epsilon.[n] is a noise term. 

Detailed Description Text (10) : 

After data acquisition, the raw data is pre-processed as ID data. That is, once 
non -water suppressed reference data is collected 114, averaging is used to yield an 
averaged set of reference data, r[n], where n is an index representing each sampled 
complex data point. This raw data averaging is done for each coil. 

Detailed Description Text (11) : 

The reference data, r[n], undergoes DC mixing 116, zero phasing 118, linear phase 
correction 120 and phase spline smoothing 122, to determine a phase correction 
vector 124, c[n]. Additionally, the reference data, r[n], is used to determine 
weighting 125 for the data of a respective coil. Simply, weighting 125 is carried 
out to determine data influence from a particular coil on a combined spectrum or 
image. Specifically, if it is determined that the data of a particular coil should 
be used, the reference data, r[n], is also used to assign weight to the receiver 
channels . 

Detailed Description Text (12) : 

The weighting 125 can be determined through multiple considerations. For a 
particular coil, a maximu m magnitude for the averaged reference data, r[n], can be 
used as a particular weighting factor, w.sub.i. Additionally, channels that have a 
significantly weak signal relative to the channel with the strongest signal are not 
used for the combined results and assigned a weight factor of zero. Specifically, 
in multi -channel MRS, regions of interest close to receive coil elements benefit 
from improved signal -to -noise ratios . Therefore, data received from regions of 
interest close to receive coil elements are appropriately weighted to take 
advantage of the improved signal -to -noise ratios . The signal characteristics from 
each coil depend on a number of factors including the orientation of the coil with 
respect to the B.sub.O field, the proximity of the coil to the volume generating 
the signal, coil loading, coil-to-coil coupling effects, and the permeability and 
permitivity of the medium through which the radio- frequency signal travels prior to 
being received by the coil elements . Accordingly, the current invention weights 
^^ta received to take advantage of the improved signal -to -noise ratios . Once 
weighted, the data is then normalized according to: ##EQU2## 

Detailed Description Text (14) : 

Water -suppressed data is also collected from each coil element and the data from 

each coil element is averaged 112 to obtain an averaged water - suppressed data set, 
s[n]. A phase correction vector 124, c[n], is then applied 127 to the reference 
data set, r[n], and the water -suppressed data set, s[n] at 126. The phase 
correction vector, c[n], may be defined by: 

Detailed Description Text (16) : 
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wherein . omega .. sub . m is the largest frequency component of the non -water 
suppressed data set, s[n], . omega .. sub . Ip is the number of phase cycles in the non- 
water suppressed data set, s[n], and .phi.. sub. s [n] is the phase spline smoothing 
factor from an unwrapped phase of: 

Detailed Description Text (17) : 

Applying the phase correction vector, c[n], removes residual water signal from the 
water -suppressed data, s[n], and also removes distortions due to eddy currents. 

Specifically, the correction vector, c[n], is multiplied by the reference data set, 
r[n], and the water -suppressed data set, s[n] as: 

Detailed Description Text (19) : 

respectively. Water subtraction 128, in which a scaled version of the reference 
data set is subtracted from the water - suppressed data set of each receiver channel, 
is performed and a Fourier Transform is applied 130 to the phase -corrected, water - 
suppressed data set, S . sub . corrected [n] , with residual water removed. Prior to 
computing a frequency spectrum, S . sub . results [k], the data is zero-padded for 
greater resolution. The results of the Fourier Transform 130, for each receiver 
coil, where the number of receive coils is L, is combined 132 by summing the 
computed results from the frequency spectrum S . sub . results [k], wherein 

Detailed Description Text (21) : 

where a. sub. scale is a constant representing the ratio of largest magnitudes of 

frequency components from the water suppressed and non -water suppressed data sets: 
##EQU3## 

Detailed Description Text (23) : 

^ summation of the independently processed receiver coil data sets is performed to 
yield a single superset of MRS results that takes into account each coil ' s 
respective weight, such that: ##EQU4## 

Detailed Description Text (24) : 

The absorption spectrum 134, S . sub . combined [k] is then plotted for the combined 
real values of the complex Fourier Transform of the phase -corrected, water- 
suppressed signal with residual water removed. As such, a ID absorption spectrum is 
displayed on a monitor or other area interface that is an assimilation of the data 
sets from each of the multiple receiver coils influenced by each coil's weighting . 
The absorption spectrum shows the chemical shift of various metabolites in parts 
per million (ppm) and can be used for the quantification of metabolites in vivo. 
The resulting ID spectral plot includes peaks from which a clinical diagnosis may 
be made without a separate interpretation of the results of each receiver coil. 

Detailed Description Text (25) : 

For ID plots representing the absorption spectrum, the peak area or the peak height 
representing a particular frequency component, can be calculated and compared to 
other peaks . For example, referring to FIG. 3, a ID plot 140 representing a 
combined absorption spectrum for multiple receiver channel data having a number of 
peaks is shown. While such a plot is not in and of itself novel, to produce such a 
single plot from multiple receive channels as described herein is novel. When 
reviewing the ID plot or spectrum 140, a medical professional often uses a creatine 
peak 142 as a reference for in vivo MRS studies. In some cases, the mere presence 
of a particular peak may indicate a clinically significant diagnosis. For example, 
lactate 144 which, in the example, shows up as a doublet, should never occur in a 
normal human brain. Therefore, the presence of lactate 144 would potentially 
indicate damage from acute stroke. As a further example, if the NAA peak 146 was 
lower and the creatine peak 142 and the choline peak 148 were higher, the plot 140 
would be indicative of a patient stricken with Alzheimer's disease. 

Detailed Description Text (26) : 

Additionally, beyond the Fourier transform, several other nonparametric techniques 
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can be used for generating results from MRS studies. Two nonparametric techniques 
that can be used for MRS analysis, which employ adaptive filter bank approaches, 
are the ID Capon method and the ID APES method. The filter bank approach to power 
spectral estimation employs a band-pass filter with a fixed bandwidth to estimate 
the power spectral density (PSD) of a given frequency component. The basic 
periodogram, as described in Eqn. 11, is a filter bank approach based on the 
standard Fourier transform. ##EQU5## 

Detailed Description Text (27) : 

Capon and APES are filter bank approaches that improve the estimate of the PSD by 
creating one data-dependent band-pass filter for each spectral point being 
estimated. The band-pass filter is created by a least-squares minimization process, 
which attempts to minimize the total output power of the filter, yet pass the 
frequency component of interest unaltered. A second least-squares process is then 
implemented to estimate the amplitude of the filtered signal. Capon and APES 
provide more accurate spectral estimates with lower sidelobes and narrower spectral 
peaks than the Fourier Transform based periodogram techniques. 

Detailed Description Text (28) : 

However, ID Capon and ID APSES methods fail to take into account, at least in a 
direct way, the damping associated with each signal component in Eqn. 1. Thus, 
since the detected amplitude will be altered by the damping, it is difficult to 
accurately estimate the amplitudes of the various signal components. To alleviate 
this problem, 2D Capon and 2D APES techniques can be used to obtain a 2D absorption 
spectrum plot. 

Detailed Description Text (29) : 

One benefit of a high-resolution 2D spectral plot derived by Capon or APES is that 
peaks that are close in frequency may be separated if the associated damping 
factors are different from one another. Another benefit of the 2D spectral plots is 
that the damping information provided by Capon or APES provides an effective means 
of estimating T . sub . 2 . sup . . multidot . , the effective decay rate of transverse 
magnetization, for several metabolites of interest. Furthermore, unlike ID methods, 
the true amplitudes are not masked by damping effects. Empirical results have shown 
that Capon analysis provides a more accurate estimate of frequency and damping 
since it provides fine resolution to resolve closely spaced peaks , while APES 
provides a more accurate estimate of the amplitude for a given frequency and 
damping factor. As will be shown, 2D absorption plots may also benefit from the 
current invention through its multiple receiver channel weighting and combination 
techniques . 

Detailed Description Text (31) : 

This constraint is based upon the extent that the sinusoid of interest passes 
undistorted through the filter, while the total output energy of the filter is 
minimized. Two-dimensional Capon and 2D APES provide an estimate of both the 
amplitude and phase for each frequency and damping factor of an input signal. The 
2D Capon and 2D APES methods provide high-resolution 2D MRS results showing both 
frequency, .omega., and damping, .sigma.. In the case of a 2D proton spectroscopy, 
the 2D Capon method and the 2D APES method, each improves upon the results of ID 
proton spectroscopy . The results have more defined peaks and more precise 
amplitude . Additionally, a damping factor representing the exponential decay of the 
associated metabolite is also provided and is useful for clinical diagnosis. 

Detailed Description Text (41): 

Referring again to FIG. 2, MRS data is again acquired 110 from multiple receive 
channels of an MR system simultaneously. After data acquisition, the raw data is 
pre-processed. That is, once MRS data is acquired 110, non -water suppressed 
reference data is collected, and averaging 114 is used to yield an averaged set of 
reference data, r[n], for each receiver coil, where n represents normalized time. 
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Detailed Description Text (42) : 

The reference data, r[n], is subjected to the DC mixing 116, zero phasing 118, 
linear phase correction 120, and phase spline smoothing 122 to compute 124 a phase 
correction vector, c[n]. Additionally, the reference data sets, r[n], are used to 
determine whether data from a particular coil should be used 125 to generate 
combined results. If it is determined that data from a particular coil should be 
used, the reference data set, r[n], is also used to assign weights 125 to each 
receiver coil . 

Detailed Description Text (43) : 

The weightings 125 are used to increase the influence of particular data sets on a 
combined superset of data. Specifically, those channels with greater signal -to- 
noise ratios are more heavily weighted . Therefore, data received from regions of 
interest close to receive coil elements are appropriately weighted to take 
advantage of the improved signal -to -noise ratios . The signal characteristics from 
each coil depend on a number of factors including the orientation of the coil with 
respect to the B.sub.O field, the proximity of the coil to the volume generating 
the signal, coil loading, coil-to-coil coupling effects, and the permeability and 
permitivity of the medium through which the radio- freguency signal travels prior to 
being received by the coil elements . 

Detailed Description Text (44) : 

A weighted least-squares estimator for 2D Capon or 2D APES may also be implemented. 
Specifically, an exponential time - weigh ting factor to adjust the least-squares 
error minimization criterion for both 2D Capon and 2D APES is incorporated in such 
a way that data samples are weighted differently based on how early or late they 
occur in the time sequence. Thus, signal elements are preferentially weighted in 
time sequentially to take advantage of improved signal -to -noise characteristics 
earlier in the signal . The exponential time weighting factor used to adjust the 
least-squares error minimization criterion in 2D Capon and 2D APES is a novel 
technique used for processing MRS data and is thus part of the invention. 

Detailed Description Text (45) : 

As previously described, water -suppressed data sets are also collected 110 from 
each coil element and the data from each coil element is averaged to obtain a 
number of averaged, water - suppressed data sets, s[n] 112. Data acquired during the 
reference acquisition is used to compute a phase correction vector 124 that is 
applied to the water -suppressed data 126. The reference data is also used to remove 
residual water. The phase correction vector, c[n] 128, is applied to the reference 
data, r[n] 127, and the water -suppressed data, s[n] 126. Applying c[n] removes 
residual water signal from the water - suppressed data sets, r[n] . Water subtraction 
128, in which a scaled version of the reference data sets is subtracted from the 
water - suppressed data sets, is performed and a nonparametric technique, 
implementing an adaptive filter bank, is used to generate the results from the 
phase -corrected, water - suppressed data sets with residual water effects removed. 
Specifically, either a 2D Capon or 2D APES technique is employed. The current 
invention performs pre-processing of raw data prior to 2D Capon or 2D APES by 
applying the previously described method of MRS data pre-processing involving phase 
correction and residual water removal. 

Detailed Description Text (47) : 

The frequency spectrum .omega, versus damping .sigma. results for each receiver 
coil, including each receive coil's respective weighting, as previously described, 
are then combined 132. If a 2D Capon technique is employed the combination results 
in a supserset of data 154, C . sub . combined [. sigma .,. omega .] , defined by: ##EQU11## 

Detailed Description Text (49) : 

From the 2D Capon method and 2D APES method, 2D absorption spectrums, 

C . sub . combined [. sigma .,. omega .] , and A. sub . combined [. sigma .,. omega .] 140, are 
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displayed, respectively, for data acquired over multiple receiver channels. These 
absorption spectra allow a reviewing medical professional to interpret the data 
received by each receiver coil simultaneously thereby increasing diagnosing time 
and accuracy. 

Detailed Description Text (50) : 

Referring now to FIG. 4, a combined 2D absorption spectrum 158 is shown. When 
compared with the ID absorption spectrum of FIG. 3, FIG. 4 shows improved peak ~to~ 
noise ratios vs. the MRS results from FIG. 3. The ratios of peak heights vs. the 
standard deviation of a noise region (near 0.0 ppm) are provided for several peaks 
of interest. This distinction becomes more evident when comparing peak -to-noise 
ratios for 2D Capon or 2D APES with each peak of a corresponding noise region with 
the same value of damping, .sigma.. Again, by compiling a single absorption 
spectrum, rather than requiring a medical professional to review an absorption 
spectrum from each receiver coil independently and individually, the burden on the 
reviewing medical professional is lowered and typically makes the review process 
considerably more efficient. The combined 2D absorption spectrum 158 allows a 
reviewing medical professional to interpret the data received by each receive coil 
simultaneously. 

Detailed Description Text (51) : 

The invention described above is a technique for multi-channel MRS whereby MR 
signals are simultaneously acquired from multiple coils and individually processed 
to generate multiple sets of MRS results. Furthermore, the multiple sets of MRS 
results are combined to form a single superset of MRS results, such that the MRS 
results are displayed as a proton MRS absorption spectrum from the single superset 
of MRS results. It is contemplated that the above architecture and technique can be 
embodied in any modality of a medical imaging scanner capable of providing data 
required for MRS. 

Detailed Description Text (52) : 

It is further contemplated that the above technique be embodied in a MR I system 
having a plurality of gradient coils positioned about a bore of a magnet to impress 
a polarizing magnetic field and an RF transceiver system and an RF switch 
controlled by a pulse module to transmit RF signals to an RF coil assembly to 
acquire MR images. The MR I system also includes a computer programmed to acquire MR 
signals from multiple receiver channels and process the MR signals from each 
channel independently. The computer then combines the processed MR signals from 
each coil into a single set of MRS results and displays a proton MRS absorption 
spectrum from the single composite set of MRS results. 

Detailed Description Text (53) : 

It is additionally contemplated that the above technique be embodied in a computer 
program including instructions which when executed by a computer cause the computer 
to process MR data from the multiple coils independently. The computer is also 
caused to combine the processed MR data into a single composite set of MRS results 
and display the set as a single MRS absorption spectrum . 

Other Reference Publication (3) : 

Stoica, P. and Sundin, T. Nonparametric NMR Spectroscopy. Journal of Magnetic 
Resonance 152, pp. 57-69, 2001. 

Other Reference Publication (4) : 

L.L. Wald et al . , Proton Spectroscopic Imaging of the Human Brain Using Phased 
Array Detectors, MRM, Oct. '95, pp. 440-445, vol. 34, No. 3. 

CLAIMS : 

1. A method of multi -channel MR spectroscopy (MRS) comprising the steps of: 
simultaneously acquiring MR signals from multiple coils ; processing the MR signals 
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individually to generate multiple sets of MRS results; combining the multiple sets 
of MRS results to form a single superset of MRS results; and generating and 
displaying a proton MRS absorption spectrum from the single superset of MRS results 
for clinical inspection. 

3. The method of claim 1 wherein the step of processing includes steps of: 
averaging individual frames of non -water suppressed reference data to obtain a non- 
water suppressed reference data set^ r, for data from each coil ; averaging 
individual frames of water suppressed data to obtain a water suppressed data set ^ s 
[n] , for data from each coil; applying a phase correction vector to each data set r 
[n] and s[n]; and removing residual water signal from the water suppressed data set 
to generate a corrected MRS signal. 

4. The method of claim 3 wherein the steps of combining includes the steps of: 
determining a maximum magnitude of the non -water suppressed data set for data from 
each coil ; from the maximum magnitude for data from each coil determining a 
weighting factor for data from each coil; applying a respective weighting factor to 
each set of MRS results based on the coil used to acquire the MR signals for the 
respective set of MRS results as weighted sets of MRS results; and summing the 
weighted sets of MRS results to form the single superset of MRS results. 

5. The method of claim 4 further comprising the step of normalizing the weighting 
factors . 

6. The method of claim 3 further comprising the step of zero-padding time domain 
data of the corrected MRS signal and wherein the steps of generating and displaying 
a proton MRS absorption spectrum includes the step of generating a frequency 
spectrum from the zero-padded time domain data. 

8. The method of claim 3 wherein the step of generating and displaying a single 
portion MRS absorption signal includes the step of applying a nonparametric filter 

bank technique to time domain data representing the corrected MRS signal to 
generate a 2D frequency spectrum versus damping. 

10. The method of claim 3 further compromising the step of determining the phase 
vector c[n] by: multiplying the non -water suppressed data set, r[n], by a scalar, 
A.sub.zp =r*[0]; determining a largest frequency component . omega .. sub . m, of the 
non -water suppressed data set, s[n]; determining a number of phase cycles, 
w.sub.lp, in the non -water suppressed data set, s[n]; determining a phase spline 
smoothing factor .phi.. sub. s [n] , from an unwrapped phase of 

r [n] e . sup . -i (w. sub .m +w.sub.lp)n 

generating a data windowing function .omega. [n]; and multiplying 

A. sub . zp.multidot . e . sup. -i (w. sub .m +w. sub. Ip) n+. o slashed. . sub . s . sup . (n) 

and . omega . [n] . 

11. The method of claim 10 further comprising the step of determining a corrected 
non -water suppressed signal, r [n] . sub . corrected =c [n] .multidot . r [n] , and 
determining a corrected water suppressed signal, s [n] . sub . corrected =s 

[n] .multidot . r [n] . 

12. The method of claim 11 wherein the step of generating and displaying a single 
proton MRS absorption spectrum includes the steps multiplying s . sub . corrected [n] 
by a scalar, a. sub. scale, and subtracting r[n] from the product wherein a. sub. scale 
includes a constant representing a ratio of largest magnitude of frequency 
components of the water suppressed and non -water suppressed data sets, sand r[n], 
respectively. 

13. A magnetic resonance imaging (MRI) system having a plurality of gradient coils 
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positioned about a bore of a magnet to impress a polarizing magnetic field and an 
RF transceiver system and an RF switch controlled by a pulse module to transmit RF 
signals to a multi- receiver channel RF coil assembly to acquire MR images; and a 
computer programmed to: acquire MR signals from multiple receiver channels; process 
the MR signals from each receiver channel independently; combine the processed MR 
signals from each receiver channel into a single composite set of MRS results; and 
display a proton MRS absorption spectrum from the single composite set of MRS 
results . 

14. The MR I system of claim 13 wherein the proton MRS absorption spectrum includes 
a single composite 2D proton MRS absorption spectrum . 

15. The MR I system of claim 13 wherein the computer is further programmed to assign 
weight factors to the MR results from each coil before the MR results are combined. 

16. The MR I system of claim 15 wherein a set of MRS results is assigned a weight 
factor of zero if the acquired MR signal is determined to be weak as compared to 
other MR signals acquired. 

17. The MR I system of claim 15 wherein the computer is further programmed to 
normalize the weight factors. 

18. The MR I system of claim 13 wherein the computer is further programmed to apply 
a nonparametric filter bank technique to time domain data representing the 
corrected MRS signal to generate a frequency spectrum versus damping. 

19. The MR I system of claim 13 wherein the computer is further programmed to: 
average individual frames of non -water suppressed reference data to obtain a non- 
water suppressed reference data set, r, for data from each coil ; average individual 
frames of water suppressed data to obtain a water suppressed data set, s, for data 
from each coil ; and apply a phase correction vector to each data set r[n] and s[n] 
to remove residual water signal from the water suppressed data set to generate a 
corrected MRS data set for data from each coil. 



20. The MR I system of claim 13 wherein the computer is further programmed to: 
determine a maximum magnitude of the non -water suppressed data set for data from 
each coil ; determine a weighting factor for each coil from the maximum magnitude 
for data from each coil; apply the weighting factor to each set of MRS results 
based on the coil used to acquire the respective set of MRS results as weighted 
sets of MRS results; and combine the weighted sets of MRS results to form the 
single superset of MRS results. 

21. A computer readable storage medium having stored thereon a computer program to 
generate a single composite MRS spectrum from MR data acquired from multiple coil 
elements that comprises a set of instructions, which when executed by a computer, 
causes the computer to: acquire MR data from multiple coil elements independently; 
combine the processed MR data into a single composite set of MRS results; and 
display a single MRS absorption spectrum as a function of the combined set of MRS 
results acquired from the multiple coil elements. 

22. The computer program of claim 21 further comprising instructions to process the 
MR data from multiple coil clients independently. 

23. The computer program of claim 22 further comprising instructions to assign 
weight factors to the MRS results from each coil before the processed MRS results 
are combined . 

24. The computer program of claim 23 having further instructions to assign a weight 
factor of zero if an acquired MR signal is determined to be weak when compared to 
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other acquired MR signals. 

25. The computer program of claim 23 having further instructions to preferentially 
weight each data set based on data acquisition time. 

27. The computer program of claim 21 having further instructions to apply a non- 
parametric filter bank technique to time domain data representing the corrected MRS 
signal to generate a frequency spectrum versus damping. 
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